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Abstract 

We study the effective quark mass induced by the finite separation of the domain walls in 
the domain-wall formulation of chiral fermion as the function of the size of the fifth dimension 
(L s ), the gauge coupling (3 and the physical volume V. We measure the mass by calculating the 
small eigenvalues of the hermitian domain-wall Dirac operator (-Hdwf^o)) in the topologically- 
nontrivial quenched SU(3) gauge configurations. We find that the induced quark mass is nearly 



(N 



independent of the physical volume, decays exponentially as a function of L s , and has a strong 
dependence on the size of quantum fluctuations controlled by (3. The effect of the choice of the 
\Q | lattice gluon action is also studied. 
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Simulating massless or near-massless fermions on a lattice is a serious challenge in numer- 
ical quantum field theory. The origin of the difficulty can be traced to the well-known no-go 
theorem first shown by Nielsen and Ninomiya, which states that one can not write down a 
local, hermitian, and chirally-symmetric lattice fermion action without the fermion doubling 
problem [|T|. Hence to have chirally-symmetric fermions on a lattice, one must use nonlocal 
actions in which the coupling between lattice sites do not identically vanish even when the 
separation becomes large. This implies that the lattice simulation with chiral fermions is 
necessarily more expensive than, e.g, the standard Wilson fermion in which only the nearest 
neighbor coupling is involved. 

One of the lattice chiral fermion formalisms that have been studied extensively in recent 
years is the domain- wall fermion, first formulated by D. Kaplan || and later modified for 
realistic lattice simulation by Shamir ||. In the domain- wall construction, one introduces 
an extra fifth dimension s with a finite extension. After discretization, the fifth direction has 
L s number of lattice sites. If we put the same four-dimensional gauge configuration on every 
four- dimensional s slices, the five-dimensional massive theory admits a four- dimensional 
effective theory in which a left-handed chiral fermion lives near the s = slice and a right- 
handed one near s = L s — 1. Integrating out the heavy modes ||, ||, ||, one obtains an 
effective four-dimensional chiral theory in the limit of L s — > oo. 

For finite L s , however, the two chiral modes can couple to produce an effective quark 
mass. Strong gauge field fluctuations can induce rather strong coupling, and hence rather 
large quark mass. This quark mass is expected to decrease exponentially as L s — > oo with 
possible power law corrections. To gain a quantitative understanding of this induced mass, 
Columbia and other groups have considered several different ways to measure it. One way 
is to study the behavior of pion mass as the function of the wall separation L s || [J. The 
problem with this is that the pion mass may not vanish as L s — > oo due to the finite volume 
effects, and hence the non-zero pion mass cannot be entirely attributed to the quark mass 
effect. Another way is to use the Gell-Mann-Oakes-Renner relation |TD|, |TTJ. However, 



this relation assumes a more complicated form in the domain-wall fermion formalism. The 
effective mass can be also measured by studying the anomalous contribution in the axial 
Ward-Takahashi(WT) identity g || |T3 



In Ref. |14[], we proposed to measure the induced effective quark mass by considering 
the eigenvalue of the hermitian domain-wall Dirac operator in the presence of the instanton 



background. In the L s — > oo limit, the lattice version of the Atiyah-Singer theorem [T5|, [L6 
guarantees the existence of exact zero modes. For finite L s , the lowest eigenvalues of the 
Dirac operator are not zero. We take the average of these would-be-zero eigenvalues as 
the effective quark mass. Our previous study for 150 gauge configurations at (3 = 6.0 
yields an effective quark mass 0.0074 ± 0.0007 at L s = 8, 0.0022 ± 0.0003 at L s = 12, and 
0.0008 ±0.00013 at L s = 16. The result is qualitatively consistent with those obtained from 



other methods IP, 12, 13 



In this paper, we report a more systematic study of the effective quark masses along 
this direction. In particular, we would like to understand the effects of different lattice 
sizes, the coupling constant (3, and the form of gluon actions. Crucial to the size of the 
effective mass is the near-zero eigenvalues of the hermitian Wilson-Dirac operator at the 
fixed domain-wall height mo. These near-zero eigenvalues require large L s to project out the 
correct effective Dirac operator. While we find that the effective mass is largely independent 
of the lattice size, it is a sensitive function of f3 and the form of gluon actions. For a strong 
coupling (small (3) the effective low-energy theory is recovered only at very large L s . We 
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find that the exponential decay rate for the effective mass is well described by the density 
of zero eigenvalues of the hermitian Wilson-Dirac operator. For an improved gluon action, 
the quantum fluctuations are strongly reduced and hence the domain-wall formalism works 
much more efficiently 

The domain- wall fermion was first introduced by D. Kaplan [0], based on an interesting 
observation that under certain conditions a five- dimensional massive fermionic theory has 
an effective four- dimensional massless fermion. For practical lattice simulations, Shamir 
considered a five-dimensional lattice having a finite extension in the fifth direction s (s = 
0, 1, • • • L s — 1). The lattice domain-wall fermion action is a five- dimensional Wilson action 
at r = 1, 



£„-l_ 

x.s 



x,s=0 



(mo - 5)ijj XtS + Pr^ x ,s-i + Pl^x,s+i 



o - ln)Ux^x+^s + (1 + lf^Ul-f^x-fi,s) 



2 , 



(1) 



where m is the negative of the conventional Wilson mass and Pr,l = |(1 ±75) is the chiral 
projection operators. The physical significance of is the "height of the domain wall". 
Like the inverse of the lattice spacing 1/a (taken to be 1), it acts as a heavy mass scale of 
the theory. The gauge fields live on the four- dimensional sub-lattices and are independent 
of s. Dirichlet boundary conditions are applied to the fermion fields at s — — 1 and L s . For 
the sake of simplicity, we omit the Pauli-Villars fields which can be introduced to cancel the 
bulk of the heavy modes of the theory. The interested reader can consult Ref. H for details. 

As shown in Ref. ||, the above construction in the absence of gauge fields has a low- 
energy effective theory with two massless chiral fermions near s = and L s — 1 in the limit 
of L s — > 00. Thus, a single-flavor massless Dirac field q(x) can be constructed from the 
domain- wall fermion field ip(x,s) as 

q(x) = P L ij(x, s = 0) + P R ij(x, s = L s - 1) . (2) 

An explicit fermion mass, rrif, for the chiral modes can be introduced through the modified 
boundary conditions, 

PliPx,l 3 = -mfP L ipx,o , Pr^x-i = -mfPn?Px,L s -i ■ (3) 

Substituting this into Eq. (0), we obtain a generic fermion mass term. 

To study the low-energy aspects of the domain-wall fermion action, we introduce the 
hermitian domain- wall Dirac operator, if DWF (m ) = -Rs75-Ddwf(^o) > where R s denotes the 
reflection in the fifth direction. Because of the low-energy property mentioned above, we 
expect that the low-lying spectrum of H DWF (m ) resembles that of a massless Dirac particle 
and contributes dominantly to low-energy physical observables. [Unlike the four- dimensional 
Wilson-Dirac operator Dw(mo) which has a simple hermitian-conjugation property, we find 
that the domain- wall Dirac operator Ddwf(wo) has a peculiar spectrum with obscure phys- 
ical significance.] 

A more direct way of studying the low-energy physics of the domain-wall fermion is to 
integrate out all heavy modes. This has been done by a number of authors [|], [5], || 0]. The 
result has a very simple form 

S cS = £g(aOD£(a 5 ,moMy) , (4) 

x,y 
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Crossing points vs. effective masses 

L =8,12,16, m eff =7.4(0.7),2.2(0.3),0.80(0.13)xl0" 3 , (3=6.0, 8 4 lattice 
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FIG. 1: Effective quark mass induced by domain walls for the Monte Carlo configurations at 
j3 = 6.0, 8 4 lattice. L s is the number of lattice sites in the fifth direction. 

where we have made the lattice spacing in the fifth direction 05 explicit and [0] 

1 + 7 5 tanh(^#(a 5 ,m ) 

Df s (a 5 ,m ) = 777777} ; 

1 - 7 5 tanh {-fH(a 5 ,m ) 

H(a 5 ,m ) = — In . 5 

a 5 \1 - a 5 P L H w (m ) J 

iJ\v( m o) — 75-Dw( m o) is the hermitian Wilson-Dirac operator. The four- dimensional effec- 
tive Dirac operator Df^(a5, mo) can be used to calculate the propagator of the chiral fermion 
modes. If an observable involves the fermion determinant, we must take into account the 
contribution from the Pauli-Villars particles as well. The effective Dirac operator in the 
fermion determinant is 



D Ls (a 5 ,m ) = - 



I+75 tanh ( y#(a 5 , m )\ 



(6) 



In the limit of a continuous fifth dimension, i.e. — > 0, H(a<^, mo) in the above expression 
becomes 

H(a 5 = 0,m ) = H w (m ) . (7) 
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Crossing points vs. effective masses 

L =8,12,16, m = 21.2(2.8),12.0(2.4),8.1(2.1)xl0~ 3 , (3=5.85, 8 4 lattice 
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FIG. 2: Effective quark mass induced by domain walls for the Monte Carlo configurations at 
j3 = 5.85, 8 4 lattice. L s is the number of lattice sites in the fifth direction. 

And if we take the limit L s — > oo as well, 

tanh (y#) e(H) = 9(H) - 6{-H) , (8) 

the Dirac operator D^a^ = 0, mo) becomes the well-known Neuberger overlap operator [|J. 
On the other hand, if we take L s — > oo first, we have 

Doota^mo) = - [1 +7 5 e(#(a 5 ,m ))] . (9) 

The difference from the Neuberger overlap operator [|J is that if (05, mo) is now a function 
of H w (m ) (see Eq. (||)). 

The above discussion leads to a simple and natural way of defining the topology of 



lattice gauge configurations |16|]. The effective domain- wall Dirac operator at L s — > 00 
(.D 00(05, mo)) has the following interesting property. If the hermitian operator H (05, mo) has 
unequal numbers of positive and negative eigenvalues, its determinant vanishes; this means 
that Doo(a5, mo) has zero-energy eigenvalues. Inspired by the Atiyah-Singer theorem in the 
continuum spacetime, one can regard the appearance of the zero eigenvalues as the signal for 
a nontrivial topology of the lattice configuration. Since the number of positive and negative 
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levels are the same for tuq < 0, the difference between the number of positive and negative 
levels at a positive m can be calculated by tracking the spectral flow of all eigenvalues 
of H(a5,m) between m = and m = itlq and locating the level crossings (16], IT7|, 



The topological index of a gauge configuration is given by the difference of the number of 
positive and negative energy levels of H(a 5 ,m ). The effective domain-wall Dirac operator 
D 00(^5, mo) has the same number of zero eigenvalues as the topological index. Rather than 
find the eigenvalues of H(a^,m), a computationally simpler method can be used for m < 2 
and 05 = 1 where the spectral flow of the lowest eigenvalues of the hermitian Wilson-Dirac 
operator H^{m) are tracked between m = and m = m < 2. The topological index thus 
produced is the same as that of the effective domain- wall Dirac operator .D 00(05, mo) @. 

For finite L s , however, the chiral modes of D l s (15,^0) corresponding to the gauge 
field topology do not have exact zero eigenvalue. The limit is approached exponentially 
as L s — ► 00. Since the low-lying eigenvalues of the effective Dirac operator -Di s (a 5 ,m ) 
must reproduce those of H DWF (m ), the latter has the same number of exponentially-small 
eigenvalues as the topological index. We define the finite-L,, induced fermion mass m c g as 
an average of these quasi- zero eigenvalues. 

In a previous publication [14[], we have calculated m e s from 150 configurations on an 
8 4 lattice at /3 = 6.0 and tuq = 1.8. Here again we show m e g for each of the zero modes 
in Fig. |I] where the horizontal axis is the Wilson mass m representing the locations of 
the level crossings in the eigenvalue flow of Hw(m), and the vertical axis is m e g. For 
those configurations with multi-level crossings, we relate m e g- and m by assuming that the 
spectrum flow is repulsive or levels do not cross, and this simplifies the presentation. The 
results for three different L s are shown with three different symbols. 

How does the effective mass change as the coupling grows stronger (smaller (3)1 To answer 
this question, we have generated 100 SU(3) lattice gauge configurations on a lattice size 8 4 , 
50 each at /3 — 5.85 and 5.7. We measure the eigenvalue flow of the hermitian Wilson- 
Dirac operator, and calculate the eigenvalues of -£/dwf(^o) corresponding to the nontrivial 
topology of the gauge configurations. For larger lattice spacing, quantum fluctuations are 
stronger, and some of these fluctuations can be misidentified as small size instantons. It 
turns out that they can induce strong couplings between the left-and right-handed chiral 
modes and are detrimental to the existence of the low-energy effective theory. Indeed, for 
the same L s , the effective masses are much larger at the smaller /3's than those at /3 — 6.0. 
For example, with L s = 16, m eff is 0.0008 at /3 = 6.0, 0.008 at (3 = 5.85, and 0.018 at 
(3 = 5.7. 

Moreover, for the same lattice size, physical volume is larger at smaller /3, and hence 
can house more instantons. As shown in Fig. ^| for (3 = 5.85 and Fig. |3| for f3 = 5.70, 
the total number of instantons in the 50 configurations is now 32 and 96, respectively. 
The level crossings happen at larger mo compared with Fig. 1. For (3 = 5.7, there are 
several crossings very close to m = 1.8. This is in contrast to what has been observed 
at (3 — 6.0, where the crossings occur mostly around m = 1.0 [[L4], The equal spacing 
between m e g at different L s is a clear signal for the exponential decay. However, there is 
a significant variation in the rate among all the crossings, as evident in the figures. (Note: 
After the analysis for the paper is completed, the recent work of the CP-PACS [|12j and RBC 
|T3| collaborations became available, which shows the signal for varying rate of exponential 
decay and/or nonvanishing effective mass in the L s — > 00 limit. This behavior is only seen at 
a much larger L s than those studied here. However, it is interesting to note that averaging 
eigenvalues with varying exponential rate can easily reproduce the large L s behavior of the 
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Crossing point vs. effective mass 

L=16,20,24, m eff =17.8(1.0),13.9(0.9),11.4(0.8)xl0" 3 , (3=5.7, 8 4 lattice 
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FIG. 3: Effective quark mass induced by domain walls for the Monte Carlo configurations at 
j3 = 5.7, 8 4 lattice. L s is the number of lattice sites in the fifth direction. 

effective mass observed in the aforementioned reference.) 

To see the volume dependence at a fixed mo and (3, we also measure the effective mass on 
a set of 50 configurations on an 8 3 x 16 lattice at (3 = 5.85. The total number of instantons 
is now 64, doubling that on the 8 4 lattice. We notice there are several crossings very close 
to m = 1.8. The average effective mass turns out to be essentially the same as that on the 
smaller lattice. A similar conclusion can be drawn by comparing m e ^ at (3 — 6.0 and V = 8 4 
with that at the same f3 and V = 16 3 x 24 in Ref. ||, although we note that the way of 
determining the effective mass there is quite different. This weak dependence on the size of 
the volume may allow us to extrapolate our quantitative results to large lattices necessary 
for realistic simulations. 

We also measure the effective mass on a set of 200 configurations generated from the 



1-loop Liischer-Weisz gauge action |19j with the tadpole improvement. (Crossings from only 



50 lattices are shown in Fig. O.) Similar studies using various improved gauge actions, 



including an RG improved action ]20[ are reported in |TT], [L2[]. The gauge coupling ((3 = 7.9) 
corresponds to a spacing of ~ 0.16 fm, similar to (3 — 5.7 of the Wilson action. The spectral 
flow and the domain-wall eigenvalues are studied with the same Wilson-Dirac operator. As 
shown in Fig. [5], the number of instantons as well as the distribution of the crossings differs 
significantly from the Wilson action at similar lattice spacing. Because of the decrease of 
small-scale quantum fluctuations, the probability of the crossing at mo > 1.2 is heavily 
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Crossing points vs. effective masses 

L =8,12,16, m rf =18.2(1.2), 10.2(1. 1),6.9(1.0)xl0" 3 , (3=5.85, 8 3 xl6 lattice 
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FIG. 4: Effective quark mass induced by domain walls for the Monte Carlo configurations at 
j3 = 5.85, 8 3 x 16 lattice. L s is the number of lattice sites in the fifth direction. 

suppressed, and the density of small eigenvalues of ifw(^o) is also much smaller (~ 2.7 x 
10~ 4 , compared to ~ 1.8 x 10~ 3 for (3 = 5.7 Wilson action). Therefore, the effective mass 
decreases faster as a function of L s . The average effective mass at L s = 16 is ~ 6 x 10~ 3 , 
compared to 18 x 10~ 3 from the Wilson action at (3 = 5.7. 

The origin of the induced chiral fermion mass is the finiteness of L s . More specifically, 
the hyperbolic-tangent function at finite L s is used to approximate the e function (Eq. (|8])). 
This approximation is good for large eigenvalues of H(a^,m ), but poor for small ones. 
For smaller /3, H(a 5 ,m Q ) has more small eigenvalues and hence the hyperbolic-tangent is 
a worse approximation of the e-function. This point is also reflected in the dependence of 
the exponential decay rate on the density of the zero eigenvalues of H^{mo) (p(0; mo)). We 
note that since the gauge fields are replicated along the fifth dimensional slices, the relevant 
length scale in the fifth dimension is the inverse of the rate of exponential decay (a) and 
by simple engineering dimensions is given qualitatively by the density of zero eigenvalues of 
/^(Ojmo) 1 / 3 . In Fig. || we have plotted 1/a as a function of p(0;m ). The rate a for each 
coupling is calculated by fitting m e Q = m exp[— aL s ] at different L s . The statistical errors 
are estimated by doing correlated fits to single-eliminated jackknife blocks. The inverse decay 
rates from all the configurations studied show an approximate linear scaling as p(0;mo) 1 ^ 3 . 
This suggests that the density of small eigenvalues is indeed the dominating factor for the 
exponential decay rate of domain-wall fermions. 
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Crossing points vs. effective masses 

L =8,12,16, m rf =15.9(0.6), 8.7(0.4),6.0(0.4)xl0" 3 , (3= 7.9, 8 3 xl6 lattice 
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FIG. 5: Effective quark mass induced by domain walls for the 50 Monte Carlo configurations 
generated by the 1-loop, tadpole improved gauge action IS] on a (3 = 7.9, 8 3 x 16 lattice. L s is the 
number of lattice sites in the fifth direction. 

Effective masses thus obtained for the different gauge coupling and L s are plotted in 
Fig. [7|. The data from Ref. || [K| are also included for comparison. For a given gauge 



coupling, different volumes and methods of measurements have little effect on the size of the 
effective mass as well as the rate of exponential decay. However, the change of gauge action 
affects the effective mass significantly. It is quite clear from the figure [7| that for a practical 
simulation of the domain-wall fermion, one either chooses a large (3 with the conventional 
Wilson action or an improved action keeping lattice spacing large. 

To summarize, we have studied the residual chiral symmetry breaking present in domain- 
wall fermion by measuring the eigenvalues of the hermitian domain-wall Dirac operator 
corresponding to the topology of the lattice gauge configurations. Individual eigenvalues 
for the topological zero modes show clear exponential behavior in L s . We regard these 
eigenvalues as the induced mass for the surface chiral modes at finite L s separation. 

For L s and (3, we see little variation of m e g as a function of the volume. This is in 
some sense expected because the coupling of the chiral modes between the opposite walls 
has little to do with the size of the four-dimensional slice. On the other hand, a strong 
dependence on (3 is observed. In particular, the effective mass is much larger at (3 = 5.85 
or 5.7 than that at (3 = 6.0. For the improved gauge action, the spurious fluctuations are 
reduced significantly and the L s needed to obtain a good chiral symmetry is reduced. Since 
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p (0;m ) vs. a 




0.05 0.1 0.15 0.2 

(p (0;m =1.8)) 1/3 



FIG. 6: Average coefficient of the exponential decreases as a function of p(0; mo). The extrapolation 
of the fit to the continuum limit for the Wilson gauge backgrounds is shown on the left. p(0; mo) 
for the Wilson and improved gauge action are from Ref. @ for 8 3 x 16 lattices. 

the additional computation needed for the improved action is negligible, using improved 
gluon actions may enable us to simulate domain-wall fermions with larger lattice spacing. 
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